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Abstract: The non-homogeneous cubic equation with three unknowns represented by the diophantine equation 
X +Y = 7 " ZW is analyzed for its patterns of non-zero distinct integral solutions. A few interesting relations 
between the solutions and special numbers are exhibited. 
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Notations: 

t m n : Polygonal number of rank n with size in 

S n : Star number of rank n 

Pr n : Pronic number of rank n 

J n : Jacobsthal lucas number of rank tl 

J n : Jacobsthal number of rank n 

CP m n : Centered Polygonal number of rank n with size m. 

G n : Gnomonic number of rank n 

Ky n : Kynea number of rank n 

I. Introduction 

The Diophantine equations offer an unlimited field for research due to their variety [1-3]. In particular, one may refer [4-14] 
for cubic equations with four unknowns. This communication concerns with yet another interesting equation 
X 3 +Y 3 =7 2 " ZW 2 representing non-homogeneous cubic with four unknowns for determining its infinitely many non- 
zero integral points. Various relations between the solutions and special polygonal numbers, centered polygonal numbers, 
Jacobsthal numbers and kynea numbers are exhibited. 

II. Method of Analysis 

The cubic equation with four unknowns to be solved for its distinct non-zero integral solution is 

X 3 +y 3 =l 2n Z W 2 (1) 

Introduction of the transformations, 
x = u + v, y = u - v, z = 2u (2) 

in (1) leads to 

U 2 +3v 2 =7 2 "w 2 (3) 

We present below different methods of solving (3) and thus, in view of (2), different patterns of solutions 
to (1) are obtained 

Pattern: 1.1 

Let w = a 2 +3b 2 (4) 
write 7 as 

7 = (2 + /V3)(2-/V3) (5) 

Using (4) and (5) in (3) and applying the method of factorization, define 

(u + iSv) = (2 + i S) 2n {a + iSb) 2 (6) 
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Since the complex number raised to any integer power is also a complex number, we write 

(2 + iS) Zn =A i +i&B i (7) 

Where A = -[(2 + ?V3) 2 " + (2-iS) 2 "] 
2 

B l =-^[(2 + iS) 2 " -(2-iS) 2 "] 
2iS 

Using (7) in (6) and equating the real and imaginary parts, we have 



u = Aj(a 2 -3b 2 )-B l (6ab)} 
v = A l (2ab) + B l (a 2 -3b 2 )\ 



Using (8) in (2), we get 

x(a,b) = A l (a 2 -3b 2 + 2ab) + B x {a 2 -3b 2 -6ab)} 

y(a,b) = A 1 (a 2 -3b 2 -2ab)-B 1 (6ab + a 2 -3b 2 )i (9) 
z(a,b) = A l (2a 2 -6b 2 )-B l (12ab) 

Thus, (4) and (9) represent the non-trivial integral solutions of (1) 
Properties: 1.2 

(i) x(T ,1) = A i (ky n -2) + B l (3J 2 „ - 6 j n + 6(-l)" - 2) 

(ii) w(2" ,1) = j 2n + 2 

(iii) z(n, n + 1) = -Aj (CP 8 „ + I6t 3 n - St 4 5 +5)-B l (24t 3 n ) 
Pattern: 2.1 

Write 7 as 7 = (10) 

4 

And (5 + iS) 2n =A 2 + ij3B 2 (11) 
Where, A 2 = -[(5 + iSf" + (5 - iSf" ] 

B 2 - -^[(5 + iS) 2n - (5 - iS) 2n 1 
2iS 

Using (4), (10), (11) in (3) and employing the method of factorization, we have 

u + iSv = 4- (A + i JlB 7 ){a 2 -3b 2 + 2iSab) 
2 2n 

Equating real and imaginary parts, we get 

u = -^[A 2 (a 2 -3b 2 )-B 2 (6ab)]\ 

2 I (12) 

v = —[A 2 (2ab) + B 2 (a 2 -3b 2 )]\ 

Thus, taking a = 2" A and b = 2" B the non-zero distinct integral solutions to (1) are given by 
x(A, B) = A 2 {A 2 - 3B 2 + 2AB) + B 2 (A 2 - 3B 2 - 6AB) 
F(A, B) = A 2 {A 2 - 3B 2 - 2AB) - B 2 {A 2 - 3B 2 + 6AB) 
Z(A, B) = A 2 (2A 2 -6B 2 )-B 2 (12AB) 
W(A,B) = 2 2 " (A 2 +3B 2 ) 

Properties:2.2 

0>(2M) = [37 2 „+l][j 2 „+2] 
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(13) 
(14) 



(15) 
(16) 



(ii) y(l, n) = A 2 [-2t 5n - 3 Pr„ + 3t 4 n + 1] - B 2 [S n + 3t 4 n + 2] 

(iii) x(n +l,7i) = A 2 [4Pr„ - 4t 4 n + 1] - B 2 [4t 4n + 4Pr„ - 1] 

Pattern: 3.1 

Introduce the linear transformations 

u = a + 3T,v = a-T 

Let w = a 2 + I2b 2 
Write 7 as 

7 ^ (4 + /-Vl2)(4-/Vl2) 
4 

and (4 + /Vl2) 2 " = (A 3 + /7l25 3 ) 
Where A 3 = | [(4 + i Vl2)) 2 " + (4 - /Vl2) 2 " ] 

5 3 =^^[(4 + /Vl2)) 2 ' ! -(4-/Vl2) 2 "] 
2/V12 

Using (14), (15) and (16) and employing the method of factorization, define 

2a + i4v2T = -L[A 3 + i4Y2B 3 )(a 2 -12b 2 + i2-J\2ab) 

Equating real and imaginary parts, we have 
a = -L-[A 3 (a 2 -12b 2 )-B 3 (24ab)]\ 

r = — [A 3 (2ab) + B 3 (a 2 -12Z? 2 ) 
Substituting (17) in (13), we get 

u — —^-[A 3 (a 2 -\2b 2 + l2ab)-B 3 (24ab-6a 2 +12b 2 )\ 

2 r i (i8) 

v = ^-[A 3 (a 2 -12Z? 2 -4ab)-BA24ab + 2a 2 -24b 2 )] 

2 2,1+1 J 

Replacing a by A2" +1 and b B2" +l , the coiTesponding integral solutions are given by 

x(A, B) = 2[A 3 (2A 2 - 245 2 + SAB) - B 3 (48A5 - 4A 2 + 48fi 2 )] 

y(A, B) = 2[A 3 (16A5) - B 3 (-8A 2 + 96B 2 )] 

z(A, B) = 4[A 3 (A 2 - 125 2 + 12A5) - B 3 (24A5 - 6A 2 + 725 2 )] 

w(A,5) = 2 2 " +2 (A 2 +125 2 ) 

Properties: 3.2 

(i) x(n,l) = 2[A 3 (CP 4 „ + 6Pr„ - 6t 4 n - 25) - B 3 (-t 10n + 45Pr„ - 45t 4 n + 48)] 

(ii) y(n + l,n) = A 3 (32Pr„ ) - 165 3 {t 24n + 8 Pr„ - %t 4n + 1) 

(iii) w(2\\) = j 4n+1 + 9J 2n+3 -2 

Note: .3.3 

Replacing (13) by U = a - 3T and V = a + T (19) 

And repeating the process as in pattern. 3 the corresponding non-zero distinct integral solutions to (1) are obtain a 

x(A, B) = 2[A 3 (2A 2 - 245 2 - 8A5) + B 3 (-4SAB + A 2 - 125 2 )] 

y(A,B) = 2[A 3 (-l6AB)-B 3 (3A 2 -36B 2 )] 

z(A, B) = 4[A 3 (A 2 - 125 2 - 12A5) - B 3 (24AB + A 2 - 125 2 )] 
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w(A,B) = 2 2 " +2 (A 2 +12B 2 ) 
Properties: 3.4 

(i) x(n,l)=A 3 (t 10n -13Pr n +13? 4n -48) + B 3 (t 6n - 95 Pr n + 95t 4 n - 24) 

(ii) y(n,l) = A 3 (-32Pr„ + 32t 4n ) + 65 3 (12/ 4 „ - 1) 

(Ui)z(n,l) = A 3 (t l0n -45Pr„ + 45r 4 „ -48)-5 3 (CP 8 „ + 95Pr n -95f 4n -49) 

Pattern:4. 1 

Instead of (15) we write 7 as 
?= (10 + i^2)(10 -»Vl2) 

16 _ 
And (10 + /V12) 2 " = (A 4 +ijY2B 4 ) 

Where A 4 = ^[(10+ iy/l2) 2n + (10 - iyll2) 2 " 

B 4 =— ^=[(10 + i4\2) 2 " -(10-/V12) 2 ' 
2/Vl2 

Using (14), (20) and (21) and equating real and imaginary parts, we have 

a = -!—[A 4 (a 2 -I2b 2 )-B 4 (24ab)]\ 
7 

(22) 

T = —[AA2ab) + B 4 (a 2 -\2b 2 )} 
2 4 " J 

Substituting (22) in (13), we get 

u = — l — [AAa 2 -12b 2 + \2ab)-BA24ab-6a 2 + 12b 2 )] 

2 4n+1 

v = — ^— [AAa 2 -\2b 2 -Aab)-BA2Aab + 2a 2 -24b 2 )] 

2 4,1+1 

To get a integer solution replacing a by 2" +1 A and b by 2" +1 B 

x(A, B) = 2[A 4 (2A 2 -24B 2 + SAB) - B 4 (4SAB - A A 2 + 4SB 2 )] 

y(A, B) = 2[A 4 (16AB) - B 4 (-SA 2 + 965 2 )] 

z(A,B) = 4[A 4 (A 2 -12B 2 + 12AB)- B 4 (24AB -6A 2 +12B 2 )] 

w(A,B) = 2 4n+2 (A 2 +12B 2 ) 

Properties: 4.2 

(i) x(2" ,1) = 4A 4 (ky n +2j n -U- 2(-l) n ) + SB 4 (ky„ - 14 j n -11 + 14(-1)" ) 

(ii) y(n,\) = A 4 (32t 3 n - 32t 4n ) - B 4 (-16? 4 „ + 192) 

(iii) z(n + l,n) = 4[A 4 (t 4n + lG n + S)-B 4 (CP l6n + CP 20n + 7St 4 „ - 8)] 

Note: 4.3 

Using (19) and repeating the process as in pattern.4, the non-zero distinct integral solutions to (1) are given by 

x(A, B) = 2[A 4 (2A 2 -24B 2 - SAB) - B 4 (4SAB + 4A 2 - 48£ 2 )] 

y(A, B) = 2[A 4 (-16AA) - B 4 (SA 2 - 96B 2 )] 

z(A,5) = 4[A 4 (A 2 -125 2 -12Afi)-5 4 (24AB + 6A 2 -725 2 )] 

w(A,B) = 2 4 " +2 (A 2 + 125 2 ) 



(20) 
(21) 
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Properties:4.4 

(i) x(n,l) = A 4 [4(t 6n -3Fr n + 2t 4n - 12)]- SB 4 [CP 4n + 10Pr„-lJi 4 „ -13] 

(ii) y(n + l,n) = -64A 4 t 3n + l6B 4 (t Sn + 8* 4 „ -1) 

(in)z(n + l,n) = -4A(CP 20n + \3t 4n - 1) + 24B 4 (t l2n + 2t 4n ) 



Pattern: 5.1 

(3) can be written as 

3v l"w + U p 

= = — ^ 0 

l"w-u v q 

Which is equivalent to the system of equations 
pu + 3vq-T'wp = 0 
qu — pv + 7"qw = 0 

Applying the cross-multiplication method, we get 
u=l n Oq 2 -p 2 ) 
v = -2*7" pq 
w = -p 2 -3q 2 



(23) 

(24) 
(25) 



Thus, the corresponding non zero distinct integral solutions to (1) are given by 
x = l"{3q 2 - p 2 - 2pq) 
y = l"{3q 2 - p 2 + 2pq) 
Z = 2*l"(3q 2 -p 2 ) 
w = -p 2 -3q 2 

Properties: 5.2 

(i) l"[x(7" ,1)] is a difference of two square 

(ii) x(7" 4) + y{T ,1) = 0(mod7) 

(Hi )6 * 7" [x(7" ,7" ) - y(7" ,7" )] is a nasty number 

(iv) z(7 n ,1" ) is a perfect square 

(v) w(2M) = -(7 2n +2) 

Pattern: 5.3 
(23) can be written as 

V l"w + U p 

7 w-u 3v q 

Repeating the process as in pattern.5, the non-zero distinct integral solutions to (1) are obtain as 
x = 7"(q 2 -3p 2 -2pq), y = T(q 2 -3p 2 +2pq) 

z = 2*T'(q 2 -3p 2 ), w = -(3p 2 +q 2 ) 

Properties: 5.4 

(i) x(n,l) = -1" (CP 6 n - Pr„ + t 4 n - 2) 

(ii) y(\,n) = l(CP 4n -t 4n -4) 

(iii) z(n, n + l) = 2* l(-2t 4n + G„ + 2) 
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III. Conclusion 

To conclude, one may search for other pattern of solutions and their corresponding properties 
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